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Abstract. We present spectrophotometric ISO imaging with the LWS and the CAM-CVF of the Serpens molecular 
cloud core. The LWS map is centred on the far infrared and submillimetre source FIRS 1/SMM 1 and its size is 
8' x 8'. The fine structure line emission in [O i] 63 /im and [C n] 157 /im is extended on the arcminute scale and can 
be successfully modelled to originate in a PDR with Go = 15 ± 10 and n(H2) in the range of (10 4 — 10 5 ) cm -3 . 
Extended emission might also be observed in the rotational line emission of H2O and high-J CO. However, lack 
of sufficient angular resolution prevents us from excluding the possibility that the emssion regions of these lines 
are point like, which could be linked to the embedded objects SMM 9/S 68 and SMM 4. Toward the Class source 
SMM 1, the LWS observations reveal, in addition to fine structure line emission, a rich spectrum of molecular 
lines, superposed onto a strong, optically thick dust continuum (Larsson et al. 2000). The sub-thermally excited 
and optically thick CO, H2O and OH lines are tracing an about 10 3 AU source with temperatures higher than 
300 K and densities above 10 6 cm -3 (M = 0.01 M©). The molecular abundances, X = N (mol) / N (R 2 ) , are X = 
(1, 0.1, 0.02, > 0.025) x 10" 4 for CO, H 2 0, OH and 13 CO, respectively. Our data are consistent with an ortho-to- 
para ratio of 3 for H2O. OH appears highly overabundant, which we tentatively ascribe to an enhanced (X-ray) 
ionisation rate in the Serpens cloud core (£ 10~ 18 s~ 1 ). We show that geometry is of concern for the correct 
interpretation of the data and based on 2D-radiative transfer modelling of the disk/torus around SMM1, which 
successfully reproduces the entire observed SED and the observed line profiles of low-to-mid- J CO isotopomers, 
we can exclude the disk to be the source of the LWS-molecular line emission. The same conclusion applies to 
models of dynamical collapse ('inside-out' infall). The 6" pixel resolution of the CAM-CVF permits us to see 
that the region of rotational H2 emission is offset from SMM 1 by 30", at position angle 340°, which is along the 
known jet flow from the Class object. This H 2 gas is extinguished by Av = 4.5mag and at a temperature of 
10 3 K, which suggests that the heating of the gas is achieved through relatively slow shocks. Although we are not 
able to establish any firm conclusion regarding the detailed nature of the shock waves, our observations of the 
molecular line emission from SMM 1 are to a limited extent explainable in terms of an admixture of J-shocks and 
of C-shocks, the latter with speeds of about (15-20) kms -1 , whereas dynamical infall is not directly revealed by 
our data. 

Key words. ISM: individual objects: Serpens cloud core, FIRS 1/SMM 1 - abundances - molecules - clouds - jets 
and outflows - Stars: formation 



1. Introduction 

The Serpens dark cloud is currently forming a dense 
cluster of low to intermediate mass stars (Strom et al. 
1976, Kaas 1999). This star forming complex is situated 
in the inner Galaxy, not very far from the direction to- 
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ward the Galactic Centre (£ n = 32°). At the distance of 
310 pc (deLara et al. 1991), the Serpens cloud core is only 
30 pc from the nominal galactic plane, i.e. well within the 
scale height of the molecular gas of about 80 pc (Dame et 
al. 1987). Not totally unexpected, IRAS-observations re- 
vealed intense and patchy far infrared (FIR) emission on 
a sloping background. Earlier observations had discovered 
discrete FIR sources of relatively low luminosity (Nordh 
et al. 1982, Harvey et al. 1984). 
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The empirical classification scheme developed by Lada 
& Wilking (1984), and later extended by Andre et al. 
(1993), is based on the Spectral Energy Distribution 
(SED) of the young stellar objects. In a previous pa- 
per (Larsson et al. 2000), we were discussing the SEDs 
of the submm-sources in the Serpens cloud core and for 
the dominating source, SMM 1 (also known as Serpens 
FIRS1; Harvey et al. 1984), we reached the conclusion 
that its SED classifies it as Class 0. The existence of a 
circumstellar disk has been announced by Brown et al. 
(2000). An immediate question is then, whether SMM 1 
shows any detectable or deducable spectroscopic evidence 
of disk accretion and/or of dynamical infall. Over a 0.2 pc 
(about 2') region toward the Serpens cloud core, Williams 
& Myers (2000) reported signs of infall. Observations to- 
ward SMM 1 of molecular line profiles (Mardones et al. 
(1997) and Gregersen et al. (1997)), have not so far been 
able to reveal a clear 'collapse signature'. This is possibly 
because of confusion with the known mass outflow activ- 
ity of the source (Rodriguez et al. 1989, Eiroa et al. 1992, 
McMullin et al. 1994, White et al. 1995, Davis et al. 1999, 
Williams & Myers 2000, Testi et al. 2000) and/or different 
velocity components in the could complex. On the other 
hand, Hogcrheijde et al. (1999) found, from continuum in- 
tcrfcromctric observations, the density profile of the source 
to be consistent with the theoretical expectation of a col- 
lapsing cloud. 

The earlier molecular line results were based on obser- 
vations of low excitation transitions, which are not par- 
ticularly (or not all) sensitive to the conditions expected 
to prevail in the deeper layers of the source. In this pa- 
per, we present spectral line data of the Serpens cloud core, 
both for a spatial map and for pointed deep integrations, 
which contain lines also of very high excitation. These are 
potentially better suited to 'penetrate' to regions which 
were previously hidden from view. Our interpretation of 
the results will especially focus on the evolution of this 
star forming complex. 

In Sect. 2, we reiterate the LWS observations and a 
summary is given for the data reductions, whereas a more 
detailed account is provided in Appendix A. The resulting 
line spectra are presented in Sect. 3. These results are dis- 
cussed at some depth in Sect. 4. We first exploit relatively 
simple analytical methods, which result in spatially aver- 
aged properties. These should be useful to limit the pa- 
rameter space for more sophisticated numerical modelling. 
This is done for the transfer of both continuum and line 
radiation, and these models lead to some valuable conclu- 
sions. In a summarising discussion, we make an attempt 
to bring these various pieces of information into a coher- 
ent physical picture. Finally in Sect. 5, we summarise our 
main conclusions from this work. 



2. Observations and data reductions 



2.1. The ISO observations 

A 5 x 5 spectrophotometric map in the far-infrared (FIR) 
of the Serpens cloud core was obtained with the Long- 
Wavelength Spectrometer (LWS; 43 - 197 /mi, R x = 140- 
330) on board the Infrared Space Observatory (ISO) on 
October 21, 1996. The ISO-project is described by Kessler 
et al. (1996) and the LWS is described by Clegg et al. 
(1996) and Swinyard et al. (1996). 

The formal map centre, viz. a = 18 h 29 m 50 s 29 and 6 = 
1° 15' 18"6, epoch J 2000, coincides to within 10" with the 
position of the sub-millimetre source SMM1. The point- 
ing accuracy in the map is determined as 1" (rms). The 
spacings between positions in the map, oriented along the 
equatorial coordinates, are 100". The size of the LWS- 
map is 8' x 8', corresponding to (0.7 x 0.7 = 0.5) £>§ 10 pc 2 , 
where D310 denotes the distance to the Serpens cloud core 
in units of 310 pc. 

At each map-point the grating of the LWS was scanned 
6 times in fast mode, oversampling the spectral resolution 
at twice the Nyquist rate. Each position was observed for 
nearly 15min. The centre position was re-observed half 
a year later on April 15, 1997, for a considerably longer 
integration time (24 spectral scans). 

In addition, imaging spectrophotometric data ob- 
tained with the Continuous Variable Filter (CVF; 5 - 
16.5 /mi, 32 x 32pxl, pixel-FOV = 6", R x >35) of ISO- 
CAM (Cesarsky et al. 1996) were also analysed. The de- 
tailed observing log is given in our previous paper (Larsson 
et al. 2000). 

2.2. The reduction of the ISO data 

2.3. LWS data 

The LWS data reductions were done with the interactive 
analysis package LIA. This pipeline processing was done 
in two ways. For the long wave detectors LW1 - LW5, 
the Relative Spectral Response Functions (RSRFs) of the 
most recent version (OLP 10) provided an overall better 
agreement than the older OLP 8 for the detector inter- 
calibration. However, the OLP 10 RSRFs also contained 
a number of strong features, which resulted in spurious 
'lines' in the reduced spectra (Appendix A). We decided 
therefore to use the RSRFs of OLP 8, which did not show 
this behaviour, but scaled to the OLP 10 absolute lev- 
els. For the short-wave detectors SW 1 - SW5, no such 
difference between OLP 8 and OLP 10 was apparent and 
consequently we used the latest version. A detailed ac- 
count for the first steps in the LWS reduction is given in 
Appendix A. 

Subsequent post-pipeline processing used the package 
ISAP. At each position, the individual LWS scans for each 
of the 10 detectors were examined, 'deglitched' and aver- 
aged. Corrections were applied to the 'fringed' spectra in 
the map. The fringing indicates that the emission is ex- 
tended and/or that point sources were not on the optical 
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Table 1. Observed H 2 fluxes near SMM 1 



CVF 7 fxm 



Transition 


Peak T 


Total 1 




/i a — 21 \T7 — 2\ 

(10 W cm ) 


/in — 20 tit — 2 \ 

(10 W cm ) 


0-0 S(2) 


2.9 ±0.9 


1.5 ±0.2 


0-0 S(3) 


5.2 ±0.6 


4.8 ±0.4 


0-0 S(4) 


8.8 ± 1.0 


3.6 ±0.3 


0-0 S(5) 


19.8 ± 1.2 


7.0 ±0.5 


0-0 S(6) 


6.2 ±0.8 


0.7 ±0.2 


0-0 S(7) 


15.6 ± 2.8 


5.8 ±0.6 


1-0 S(l) u 


2.5 


1.0 



Notes to the table: 

t Peak emission into a 36 arcsec 2 pixel. 

* Total emission over a source of 350 arcsec 2 . 

u Estimated from Fig. 3 of Eiroa and Casali (1989). 



axis of the LWS, i.e. the radial distance from the optical 
axis was typically larger than about 25". 

The absolute flux calibration is better than about 30% 
(Swinyard et al. 1996). Overlapping spectral regions of 
adjacent detectors were generally within 10% ('detector 
stitching' uncertainty). At 60 /im and 100 /im, the corre- 
spondance with broad-band IRAS data is even better than 
10% (Larsson et al. 2000). 

2.4. ISO-CAM-CVF data 

The ISO-CAM data were reduced with the CIA pro- 
grams using OLP8.4. This included dark current correc- 
tion, transient correction and deglitching and flat field 
corrections. The calibrations were based on OLP 5.4. The 
proper alignment of the CVF fields was achieved with the 
aid of point sources identified in both CVF and broad- 
band (7/Ltm and 14 /im) CAM observations (Fig.|l|). 



3. Results 
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Fig. 1. The from CVF frames synthesised images at 7 /im 
and 14 /im. From comparison with CAM broad-band ob- 
servations at these wavelengths and deep K-images the 
point sources were identified (see: Kaas 1999). These are 
shown by plus-signs (white and black) and were used to 
align the CVF frames with the equatorial coordinate sys- 
tem in the sky. Major SMM sources in the field are shown 
by yellow crosses. The LWS map pointings are indicated 
by the contours of the FWHM of the beam (white circles). 



3.1. Spectral imaging in H2 

Lines from purely rotationally excited H 2 are clearly dis- 
cernible in the CVF-image. In particular, all lines which 
fall in the observed spectral range, were detected and, in 
Table [j], their observed fluxes are listed. These refer to the 
3 para-H 2 lines 0-0 S(2) 12.3, S(4) 8.0 and S(6) 6.1 /mi and 
the 3 ortho-H 2 lines 0-0 S(3) 9.7, S(5) 6.9 and S(7) 5.5 /im. 
Fig. U shows maps near SMM 1 in these six lines, delineat- 
ing the spatial distribution of the rotationally excited H 2 
gas along the known jet-flow from SMM 1 (Eiroa & Casali 
1989, Hodapp 1999). In addition, somewhat weaker H 2 
emission is also observed near SMM 3, the discussion of 
which will be postponed to a future paper. 



3.2. LWS-maps of [Oi]63fim, [Cn]157nm, 
ortho-H 2 179.5 fim and CO 186p,m 

Within the 8' x 8' map observed with the LWS, the spa- 
tial distribution of the fine structure lines [O 1] 63 /im and 
[C 11] 157 /im is shown in Fig. || together with that of the ro- 
tational lines of H 2 (2i 2 — loi) (connecting to the ground 
state) and CO ( J = 14 — 13) (lowest J— transition admit- 
ted and close in wavelength to the H 2 line). 

Like the emission in [C 11] 157 /im, [0 1] 63 /im emission 
is observed in each position of our map, but with maxi- 
mum intensity at the position of SMM 1. Secondary max- 
ima are also seen along a ridge toward SMM 4/3/8 in the 
southeast-to-east. In addition, the emission is extended 
toward SMM9/S68NW. The overall [Ol] 63 /im distribu- 
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Fig. 2. CVF images toward SMM 1 in 3 ortho-H2 (upper frames) and 3 para-H-2 (lower row) lines. Offsets are in arcsec 
and relative to the centre position (0, 0) of the LWS map. The white arcs outline the FWHM contour of the LWS 
beam (70") and the white spot designates the position of SMM 1. The actual CVF pixels are square 6". 



tion is strikingly similar to that seen in the mm-regime 
in a variety of high density tracing molecules (McMullin 
et al. 2000). Very interesting is also the fact that, just in 
the northwest corner of the map, the shocked flow HH 460 
(Davis et al. 1999, Ziener & Eisloffel 1999) seems to be 
discernible. 

In [C n] 157 /im, maximum emission is seen toward S 68 
(Sharpless 1959) and the cluster in the southeast, whereas 
SMM 1 is inconspicuous in this map. To some extent, 
the [C il] 157 /im-map is an inverted image of that seen 
in [O i] 63 /im. 

The molecular emision in the 179.5 /im line of ortho- 
H2O and in the 186 /^m line of CO is strongest to- 
ward SMM 1, but displays also some extensions along the 
[0 1] 63 fim ridge, i.e. toward the northwest and the south- 
east. The spatial resolution is not sufficient, however, to 
decide whether this emission is extended or due to the 
numerous SMM-sources in the region (cf. Fig.||). 

3.3. The line spectrum of SMM 1 from 43 to 197 fim 

The continuum subtracted line spectrum of SMM 1 from 
43 to 197 /im is displayed in Fig.||, revealing a plethora of 
emission lines from both atomic and molecular species. In 
the figure, also the wavelength coverage of the ten individ- 
ual LWS detectors is indicated (cf. Sect. 2.2). The identifi- 
cation of the lines and their parameters, as obtained from 



profile fitting, are listed in Tables|| to ||. Listed are the 
identified species, the upper and lower states and the rest 
wavelength of the transition. This is followed by the ob- 
served wavelength with the measurement error and the 
difference between the observed and the rest wavelength. 
Then the line width and the line flux with fitting-errors, 
respectively, are tabulated and the error estimate for the 
line flux, i.e. AF = ^/af + a^ where <j\ is the fitting error 
and <7 C is the rms-level of the surrounding continuum inte- 
grated over one spectral resolution element. In the last two 
columns, the type of the Gauss fitting (single or multiple 
component) and the LWS detector are given. 

Out of a total of 8 measured CO lines, clear detec- 
tions are found for 7 transitions. In addition, we identfy 
5 (possibly 8) lines of ortho-H^O, and 1 (possibly 2) of 
para-H20. Further, 4 lines of OH are clearly present in 
the spectrum, whereas for two more lines this status is 
less clear. Aside from the line emission from molecules, 
lines from O and C + are also present. Possible implica- 
tions from these results will be discussed in the following 
sections. 



4. Discussion 

4.1. The atomic line spectrum: CNO 
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Fig. 3. Maps of integrated line fluxes of [Oi] 63 /im, [Cn] 157 ^m, H 2 (2 i2 - lot) and CO (14 - 13), with individual 
flux levels indicated by the bars next to each frame. The crosses refer to the positions toward which the LWS spectra 
have been obtained. 



4.1.1. [0 m] 53 fim, [0 in] 88 fim and [N n] 122 

As is evident from Table ^, only upper limits were obtained 
for the high ionisation lines [O in] 53 fim, [O in] 88 fim and 
[Nil] 122 fim. This is consistent with the luminosities de- 
rived by Larsson et al. (2000), indicating the presence of 
stellar sources generating at best only gentle UV-fields. 

We can also exclude the presence of extended strongly 
shocked regions. For instance, if associated with the fast 
moving objects of Rodriguez et al. (1989), our data imply 
that the physical scales of these shocks would be small, 
<C 1" (e.g. for v s ~200kms~ 1 and no>10 5 cm~ 3 ; see: 
Cameron & Liseau 1990; Liseau et al. 1996a). 

We can conclude that the degree of ionisation of the 
atomic gas is generally low. Lines of low-ionisation species 
will be discussed in the next sections. 



4.1.2. Extended emission: [O^lST/urn and [Oi^S/urn 

The spatial distribution of the [C n] 157 /im emission is 
shown in Fig.|^ from which it is evident that the emission 
varies within a factor of about two. The S 68 nebulosity is 



pronounced in the [C n] 157 fim line, making it likely that 
its origin is from a photondominated region (PDR), close 
to the cloud surface. 

This idea can be tested quantitatively by invoking also 
the observed [O i] 63 fim emission. The LWS subtends a 
solid angle Olws = 9.0 x 10~ 8 sr. Disregarding for a mo- 
ment the peak emission (see below), we find for unit beam 
filling the line intensities /157 = 2 x 10 _5 erg cm~ 2 s _1 
sr _1 and ^63 = 1 x 10~ 5 erg cm~ 2 s _1 sr _1 , respectively. 
Hence, the line ratio [0 1] 63 /j,m/[C 11] 157 /im is about 0.5 
(Fig. |). These data are consistent with PDR-cmission, 
where an interstellar radiation field about 10 times as in- 
tense as that of the solar neighbourhood, i.e. Go ~ 10, 
is impinging on the outer layers of a cloud with densi- 
ties in the range n(H) = (0.1 — 1) x 10 5 cm~ 3 (cf. Figs. 4 
and 5 of Liseau et al. 1999). This estimate of the strength 
of the UV field is in reasonable agreement with the FIR- 
background measured by IRAS and ISO (Larsson et al. 
2000), which would imply Go = 5 — 25. In the advocated 
PDR model, the [0 1] 145 /im line is fainter by two or- 
ders of magnitude than the [0 1] 63 ^m line. The observed 



6 Bengt Larsson et al.: The ISO-LWS map of the Serpens cloud core 




LW1 LW2 LW3 , LW4 LW5 



i i i i i i i i i i i i i i i i i i i i i i i 

100 120 140 160 180 

Wavelength [yum] 

Fig. 4. The LWS spectrum (43 to 197 /mi) of Serpens SMM 1 with line positions indicated. For clarity the strong 
continuum (Larsson et al. 2000) has been subtracted. The originally four times oversampled data have been rebinncd 
to 0.1 /im for the short wavelength (SW) and to 0.2 /im for the long wavelength (LW) spectral range. Also shown, 
below the spectrum, is the extent of the individual LWS detectors. 
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Fig. 5. Emission line intensity ratio map for 
[Oi] 63/im/[Cn] 157/im. Observed positions are identi- 
fied by crosses. The observed flux distributions of the 
individual lines are shown in Fig.pl 



line ratio, [Oi] 63/rni/[Oi] 145 /im> 12, is clearly consis- 
tent with this prediction. Finally, no detectable emission 
from higher ionisation stages would be expected. The PDR 
model would offer therefore a satisfactory explanation for 
the observed fine structure line distribution over the map. 



If this PDR emission is treated as a large scale back- 
ground and subtracted from the maps, the resulting line 
ratio toward the peak (SMM 1) increases dramatically, viz. 
to [O i] 63 /im/[C n] 157 /im = 15 ± 5. Such large ratios are 
generally not predicted by PDR models but are a com- 
mon feature of J-shocks (Hollenbach & McKee 1989). The 
residual [O i] 63 /im flux corresponds to an observed inten- 
sity Iqs = 1 x 10 _4 erg cm -2 s _1 sr _1 , more than two 
orders of magnitude below that of J-shock models (see 
Sect. 4.1.3). Interpreted as a beam filling effect, this would 
imply the size of the shocked regions to be about 4" to 5". 
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Ion id 


Transition 


A 


Wavelength 

w 

Aobs AA 




FluxxlO 19 

(W cm- 2 ) 
F AF 


b ingle/ 
Multi 
nt 


T WTQ 
L/Wb 

Detector 


[CII] 


2 ?3/2 - 2 Pi/2 


157.74 


157.72 ± 0.02 


—0.02 


0.46 ± 0.03 


2.07 ± 0.21 


0.35 


s 


LW 4 






157.72 ± 0.04 


-0.02 


0.60 Fix 


2.31 ± 0.20 


0.34 


s 


LW 4 


[01] 


3 Po - 3 Pi 


145.53 


145.53 Fix 


+0.00 


0.60 Fix 


0.37 ± 0.84 


0.90 


M 


LW 4 








145.53 Fix 


+0.00 


0.60 Fix 


0.76 ± 0.83 


0.89 


M 


LW 3 




'Pi - 3 P 2 


63.18 


63.19 ± 0.01 


+0.01 


0.27 ± 0.01 


9.5 ± 0.7 


0.74 


S 


SW 3 








63.19 ± 0.01 


+0.01 


0.29 Fix 


9.8 ± 0.4 


0.46 


S 


SW 3 








63.17 ± 0.01 


-0.01 


0.28 ± 0.01 


10.3 ± 0.9 


0.96 


s 


SW 2 








63.17 ± 0.01 


-0.01 


0.29 Fix 


10.5 ± 0.5 


0.60 


s 


SW 2 


[Nil] 


3 P 2 - 3 P1 


121.90 










0.29 


s 


LW 2 


[OIII] 


3 Pi - 3 Po 


88.36 










0.47 


s 


LW 1 
















0.32 


s 


SW 5 




3 P 2 - 3 P1 


51.82 










0.73 


s 


SW 2 



4.1.3. The [Oi]63^m emission toward HH460 

Toward the interstellar shock, HH460, the [Oi] 63 (im flux 
is not conspicuously larger than that of [Cn] 157/zm, as 
one might naively expect for shock excitation, and we 
cannot exclude the possibility that the spatial coinci- 
dence with the [O i] 63 fim emission spot is merely acci- 
dental. However, pursuing the shock idea we find that, for 
the previously inferred cloud densities, >10 5 cm~ 3 , the 
[O i] 63 /im intensity is roughly constant with the shock 
speed (a few times 10~ 2 erg cm -2 s -1 sr" 1 ; Hollenbach 
& McKee 1989). These J-shock models do also predict 
that the accompanying [O i] 145 ^m emission would not be 
detectable in our observations and that any [C n] 157 ^m 
contribution would be totally insignificant. 

The observed line intensity is 1 x 10 _5 erg cm -2 s -1 
sr -1 which, if due to shock excitation, would indicate 
that the source fills merely a tiny fraction of the LWS- 
beam (beam dilution of 2.5 x 10~ 4 ). A size of about 
1" for the [O i] 63 ^m emitting regions of the HH object 
would thus be implied, which is comparable to the dimen- 
sion of the dominating, point-like, optical knot HH460 A. 
From the observed line flux, a current mass loss rate 
from the HH-exciting source of Mi OS s = 3 x 10~ 7 Mq yr -1 
would be indicated (Hollenbach 1985, Liseau et al. 1997), 
which is at the 3% level of the mass accretion rate in the 
Serpens cloud core (Sect. 4.4). 

Based on the Lqz — Lboi calibration by Liseau et al. 
(1997), one would predict that the luminosity of the cen- 
tral source is slightly less than 0.5 Lq. No detailed in- 
formation about the exciting source of HH460 is avail- 
able, though. Based entirely on morphological arguments, 
Ziener & Eisloffel (1999) associate HH460 with SMM1, 
and Davis et al. (1999) cither with SMM 1 or with 
SMM9/S68N. The inferred luminosities of these objects, 
71 Lq and 16 L@, respectively (Larsson et al. 2000), are 
however much larger than that inferred for the putative 
source driving HH460. Evidently, the present status re- 
garding the identification of the driving source of HH 460 



is inconclusive. Proper motion and radial velocity data 
would be helpful in this context. 

4.2. H 2 , CO, H 2 and OH toward SMM1 

We can directly dismiss the PDR of Sect. 4.1.2 as respon- 
sible for the molecular line emission observed with the 
LWS, since gas densities and kinetic temperatures are far 
too low for any significant excitation of these transitions. 
Shock excitation would be an obvious option. In the fol- 
lowing, we will examine the line spectra of H2, CO, H2O 
and OH. 

4.2.1. Rotation diagram: H2 and CO 

The analytical technique known as 'rotation diagram' 
analysis is relatively simple and easy to apply to wave- 
length integrated molecular rotational line data. The ad- 
vantages and the shortcomings of this analysis tool have 
been thoroughly discussed by Goldsmith & Langer (1999). 

Assuming the lines to be optically thin and to be 
formed in Local Thermodynamic Equilibrium (LTE), one 
can derive the equation of a straight line for the molecular 
column density as a function of the upper level energy in 
temperature units. The slope of this line is the reciprocal 
excitation temperature of the levels (which in LTE is the 
same for all levels and equals the kinetic gas temperature), 
viz. 



In I, '\ k^X\ ln fN mol 



h vq g u A u i fib 



\Q(T ex )J k 



(1) 



where the symbols have their usual meaning. The left hand 
side of Eq. (1) entails the column density of the molecules 
in the upper levels. For H2, the upper level energies, 
E u , were obtained from Abgrall & Roueff (1989) and the 
Einstein transition probabilities, yl u i, were adopted from 
Wolniewicz et al. (1998). For CO, these data were taken 
from Chandra et al. (1996). The statistical weights of the 
upper levels are given by g u — (21 + 1)(2 J u + 1), where 
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Table 3. CO line measurements in the spectrum of SMM 1 



Transition Wavelength Width FluxxlO 19 ErrxlO 19 Single/ LWS 

(jjxd) (Atm) (Wcm~ 2 ) (W cm ) Multi Detector 

(J -t- J - 1) A A obs AA a x F AF fit 



14 - 


13 


186.00 


185.88 ± 0.03 


-0. 


.12 





.76 ± 0.03 


2.64 ± 0.30 


0.37 


S 


LW 


5 








185.90 ± 0.03 


-0. 


.10 




0.60 Fix 


2.32 ± 0.16 


0.27 


S 


LW 


5 


15 - 


14 


173.63 


173.60 ± 0.03 


-0. 


.03 





.81 ± 0.03 


3.34 ± 0.28 


0.35 


S 


LW 


5 








173.59 ± 0.03 


-0. 


.04 




0.60 Fix 


2.81 ± 0.21 


0.30 


s 


LW 


5 








173.58 ± 0.03 


-0 


.05 




0.60 Fix 


2.77 ± 0.50 


0.55 


M 


LW 


5 


16 - 


15 


162.81 


162.90 ± 0.02 


+0. 


.09 





.98 ± 0.02 


4.58 ± 0.19 


0.28 


S 


LW 


5 








162.92 ± 0.04 


+0. 


.12 




0.60 Fix 


3.38 ± 0.34 


0.40 


S 


LW 


5 








162.81 Fix 








0.60 Fix 


2.67 ± 0.50 


0.54 


M 


LW 


5 








162.93 ± 0.02 


+0. 


.12 





.59 ± 0.02 


3.19 ± 0.18 


0.36 


S 


LW 


4 








162.93 ± 0.04 


+0. 


.12 




0.60 Fix 


3.21 ± 0.10 


0.33 


S 


LW 


4 








162.81 Fix 


+0. 


.00 




0.60 Fix 


2.86 ± 0.34 


0.46 


M 


LW 


4 


17- 


16 


153.27 


153.28 ± 0.02 


+0. 


.01 





.64 ± 0.02 


2.81 ± 0.20 


0.34 


S 


LW 


4 








153.27 ± 0.02 


+0. 


.00 




0.60 Fix 


2.71 ± 0.11 


0.30 


S 


LW 


4 


18 - 


17 


144.78 


144.73 ± 0.02 


-0. 


.05 





.57 ± 0.03 


2.55 ± 0.31 


0.44 


S 


LW 


4 








144.75 ± 0.03 


-0 


.03 




0.60 Fix 


2.63 ± 0.16 


0.36 


S 


LW 


4 








144.78 Fix 








0.60 Fix 


2.50 ± 0.22 


0.39 


M 


LW 


4 








144.78 ± 0.02 


+0. 


.00 





.50 ± 0.02 


2.67 ± 0.28 


0.42 


S 


LW 


3 








144.79 ± 0.03 


+0. 


.01 




0.60 Fix 


2.94 ± 0.23 


0.38 


S 


LW 


3 








144.78 Fix 








0.60 Fix 


2.94 ± 0.32 


0.44 


M 


LW 


3 


19 - 


18 


137.20 


137.19 ± 0.03 


-0. 


.01 





.80 ± 0.04 


3.39 ± 0.41 


0.53 


S 


LW 


3 








137.21 ± 0.04 


+0. 


.01 




0.60 Fix 


2.85 ± 0.28 


0.43 


S 


LW 


3 








137.20 Fix 








0.60 Fix 


2.83 ± 0.30 


0.45 


M 


LW 


3 


20 - 


19 


130.37 


130.35 ± 0.03 


-0. 


.02 





.68 ± 0.03 


2.48 ± 0.32 


0.48 


S 


LW 


3 








130.34 ± 0.03 


-0 


.03 




0.60 Fix 


2.33 ± 0.18 


0.39 


S 


LW 


3 


21 - 


20 


124.19 


124.21 ± 0.05 


+0. 


.02 





.48 ± 0.06 


0.98 ± 0.32 


0.48 


S 


LW 


3 








124.20 ± 0.08 


+0. 


.01 




0.60 Fix 


1.08 ± 0.23 


0.42 


S 


LW 


3 








124.10 ± 0.05 


-0. 


.09 





.52 ± 0.06 


1.17 ± 0.31 


0.42 


S 


LW 


2 








124.10 ± 0.06 


-0. 


.09 




0.60 Fix 


1.25 ± 0.18 


0.34 


S 


LW 


2 


22 - 


21 


118.58 














0.29 


S 


LW 


2 



Notes to the tabl e: 

AF = af + g% where o\ is the fitting error and a c is the rms-level of the surrounding continuum integrated over one spectral 
resolution element, i.e. 0.29 /jm and 0.60 fim for the SW and LW detectors, respectively. S = Single line component fit; M = 
multi component line fit. Underlined flux values were used in the rotation diagram analysis and are shown in Fig. [A 



I is the quantum number of the nuclear spin. Further, for 
the evaluation of the approximate partition function 



Q(T ex ) 



HcBq 



(2) 



we used the rotational constant, Bq = 59.33451cm 1 , for 
H 2 from Bragg et al. (1982). For CO, B = 1.9225 cm- 1 , 
was obtained from the data by Lovas et al. (1979). Graphs 
of Eq. (1) are shown in Figs.| and [| fitted to the CAM- 
CVF and LWS data, respectively. 

To obtain a consistent result, the H2 data need to be 
corrected for the foreground extinction. Using the data of 
Ossenkopf & Hcnning (1994; model for thin ice coating, 
n — 10 5 cm~ 3 , t = 10 5 yr), an extinction correction of 
Ay = 4.5 mag resulted in a total column density of warm 
H 2 gas of N(H 2 ) = (9.5 ± 2.0) x 10 18 cm~ 2 . The rotation 
temperature is T ex = (1060 ± 100) K and an ortho-to-para 
ratio (nuclear spin state population) of o/p = 3 is implied 
by these data. 



In Fig. [7], ground-based CO data from the literature 
were added for lower lying transitions. Evidently, the high- 
J distribution appears markedly different from that of 
the low- J lines. If these latter lines were truly optically 
thin, they could originate in extended cloud gas, where 
T low _j = (62 ± 6) K, of column density iVi ow _j = (2.0 ± 
0.2) x 10 17 cm -2 . Seemingly in contrast, the LWS data 
identify gas at a characteristic temperature of Ihigh-,7 = 
(320 ± 20) K, with an LTE-column density of N hi gh-J = 
(1.9 ±0.2) x 10 15 cm- 2 . 

These results are based on ad hoc assumptions, i.e. 
that of unit beam filling and of low optical depth in the 
lines, potentially underestimating the column densities, 
and that the level populations are distributed according 
to their LTE values. LTE may be a reasonably good as- 
sumption for the low- J lines. Regarding CO, it is however 
questionable to what extent these are optically thin. On 
the other hand, low opacity may come close to the truth 
for the high-J lines, but LTE is not at all guaranteed a 
priori for these transitions. Obviously, one needs to check 
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Transition Wavelength Width FluxxlO 19 ErrxlO 19 Single/ LWS 

(/im) (/xm) (W cm" 2 ) (W cm" 2 ) Multi Detector 







A 


Aobs 


AA 




0"A 






F 






AF 


fit 






ortho-H 2 






























2 2 i - 


2l2 


180.49 


1 on /i n i_ a no 

180.49 ± 0.1)3 


+0.00 





A d _l_ A AO 

45 ± O.Oo 





53 


± 





11 


0.24 


b 


t xx r 
LW 


5 








18U.49 ± U.U5 


1 n An 
+0.00 




n en 1 / ; . . 
U.bU r IX 


n 
U 


OU 


± 


n 



A 1 

Ul 


n 00 
0.23 


O 
O 


LW 





2l2 - 


loi 


1 7n co 
1 <9.5o 


1 /9.55 ± U.U4 


1 a no 
+0.02 


n 



0( ± 0.04 


t 


1 A 

40 


1 

± 


n 



24 


n 00 
0.32 


O 
O 


t \mT 











179.55 ± U.U4 


+0.02 




0.60 Fix 


1 


43 


± 





13 


0.25 


c 
O 


T XX T 
LW 


5 


3o3 — 


2l2 


1 *7 A CO 
1/4. DO 


1 ta czo _i_ a no 
1 f4.0o ± U.Uo 


1 a aa 
+0.00 


n 



ta _i_ n no 
(0 ± O.Oo 



2 


1 r 

15 


1 

± 


n 



1 

24 


n 00 
0.33 


c 



T WT 
LW 


r 










1 ta go _i_ a no 
1 ( 4.UO ± U.Uo 


1 n nn 
+0.00 




n en 1 / ; . . 
U.bU r IX 


1 
t 


9( 


1 

± 


n 



1 4 

14 


n ofi 
0.26 


c 



T XX?" 
LW 


r 










174.63 ± 0.06 


+0.00 




O.bO 1 IX 


1 


93 


± 





66 


0.71 


1\1 


T X X T 

LW 


5 


330 — 


3 2 i 


1 a n 
136.49 


1 oc co _i_ n An 
13b. 52 ± U.09 


1 n no 
+0.03 


n 



c A 1 n 1 1 
b4 ± U.ll 


n 



r 

81.) 


1 

± 


n 



1 

34 


n a 
0.48 


O 
O 


T XX7" 

LW 


• > 

3 








13b. 51 ± 0.0b 


1 n no 
+0.02 




U.uO r IX 


n 



1 
81 


± 


n 



11 


a o c 
0.35 


c 



T \TiT 
LW 



3 








IOC CI T7^,. 

lob. i3l rrx 


1 n no 
+0.02 




n ca d;.. 
U.uO r IX 


n 



( ( 


1 

± 


n 



82 


n 00 
0.88 


M 


T XX?" 
LW 



3 


4l4 - 


3o3 


113.54 


i i o c _i_ n n /i 
113.54 ± U.04 


1 n nn 
+0.00 


n 



c 1 _l n n 
61 ± 0.04 



2 


AC 

9b 


1 

± 


n 



O A 

39 


n co 
0.52 


c 



LW 



2 








1 n w + (i n 1 ! 

-L-LO.O^i m U.UO 


-i-o nn 

Tw.UU 




n fin Fix 

U.UU -T IX 


2 


<)3 

.70 


-4- 


n 
u 


21 


0.41 


G 




T W 


2 


221 - 


lio 


108.07 


108.13 ± 0.05 


+0.06 





39 ± 0.02 


1 


29 


± 





40 


0.50 


S 


LW 


2 








108.18 ± 0.07 


+0.11 




0.60 Fix 


1 


58 


± 





31 


0.45 


S 


LW 


2 








107.83 ± 0.03 


-0.25 





51 ± 0.05 


1 


60 


± 





24 


0.52 


S 


LW 


1 








107.84 ± 0.05 


-0.24 




0.60 Fix 


1 


79 


± 





22 


0.51 


S 


LW 


1 


321 - 


2l2 


75.38 


75.40 ± 0.03 


+0.02 





29 ± 0.03 


2 


05 


± 





51 


0.60 


S 


sw 


4 








7CT A A _l_ A AO 

75.40 ± O.Oo 


+0.02 




0.29 Fix 


2 


03 


± 





28 


0.42 




nil 7 

sw 


4 


4 4 i - 


330 


49.34 


49.32 ± 0.02 


-0.02 





25 ± 0.02 


1 


69 


± 





29 


0.37 


s 


SW 


1 








49.32 ± 0.02 


-0.02 




0.29 Fix 


1 


81 


± 





18 


0.29 


M 


sw 


1 


para- 


H 2 






























220 — 


111 


100.98 


100.66 ± 0.05 


-0.32 





57 ± 0.04 


2 


53 


± 





68 


0.82 


s 


LW 


1 








100.66 ± 0.06 


-0.32 




0.60 Fix 


2 


60 


± 





37 


0.58 


s 


LW 


1 


322 — 


2u 


89.99 


89.99 ± 0.05 


+0.00 





44 ± 0.04 


2 


71 


± 





76 


0.83 


s 


SW 


4 








90.01 ± 0.04 


+0.02 




0.29 Fix 


2 


02 


± 





42 


0.53 


s 


SW 


4 




2000 4000 6000 8000 

E„ P [K] 



Fig. 6. Rotation diagram for the ortho-H2 (filled sym- 
bols) and para-H2 (open symbols and upper limit) lines 
observed toward the flow from SMM 1 . Linear regression 
fits to extinction corrected data is shown by the full drawn 
line. The physical parameters with their formal errors are 
given in the figure (see also the text). 

how well these assumptions are justified. In the next sec- 
tions, this will be addressed by employing first a method 
based on the Sobolev approximation and then a full Monte 



Carlo calculation, including gradients for both density and 
temperature. The latter method takes any (previously ne- 
glected) beam dilution effects directly into account. 

4.2.2. Large Velocity Gradient models: CO 

In the Large Velocity Gradient model (LVG) opacity ef- 
fects in the lines are explicitly taken into account by intro- 
ducing the photon escape probability formalism. Crudely 
speaking, the critical density of the transition, n cr j t = 
A U \/C U \(T), can be lowered by means of an effective 
Einstein-probability, A u \f3 csc , where (3 CSC is in the range 
to 1 for infinite and zero optical depth, respectively.^] 
This can effectively 'delay' line saturation. For illustrating 
purposes, (3 CSC ~ l/ri ine but, in general, (3 CSC is geometry 
dependent (Castor 1970). 

For under-resolved sources, an ambiguity can arise 
from the fact that hot and tenuous models may be indis- 
tinguishable from cool and dense ones. However, assuming 
that the rotation diagram analysis can provide an estimate 
of the kinetic gas temperature, LVG can be used to deter- 
mine the average density of the emitting region. This is 
shown in Fig.||, for the resulting logn(H2) = 6.2 ± 0.2 

The collision rate coefficients for CO, Cy, have been re- 
evaluated and extended to higher J values in Appendix B. 
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Table 5. OH line measurements in the spectrum of SMM 1 





2 II Transition 




Wavelen 


gth 




FluxxlO 19 


ErrxlO 19 


Single/ 


LWS 










w 






(W cm z ) 


(W cm z ) 


Multi 


Detector 


Q 


-n' 


J- J' 


X 


Aobs 


AA 




F 


AF 


fit 




72 


- fa 


fa - fa 


163.26 


163.26 Fix 


+0.00 


0.60 Fix 


1.39 ± 0.90 


0.93 


s 


LW 5 










163.26 Fix 


+0.00 


0.60 Fix 


1.00 ± 0.85 


0.91 


s 


LW 4 


% 

1 - 


3/ 
- /2 


5/ 3/ 
fa - fa 


119.34 


119.38 ± 0.20 


+0.04 


1.31 ± 0.11 


2.07 ± 0.49 


0.56 


s 


LW 2 










119.42 ± 0.07 


+0.08 


0.60 Fix 


1.28 ± 0.22 


0.36 


s 


LW 2 


z fa 


- /2 


7/ 5/ 
fa~ fa 


84.51 


84.43 ± 0.03 


-0.08 


0.39 ± 0.03 


2.94 ± 0.64 


0.79 


s 


LW 1 










84.44 ± 0.08 


-0.07 


0.60 Fix 


3.49 ± 0.73 


0.87 


s 


LW 1 










84.51 ± 0.02 


+0.00 


0.41 ± 0.02 


4.02 ± 0.50 


0.58 


s 


SW 5 










84.50 ± 0.02 


-0.01 


0.29 Fix 


3.26 ± 0.31 


0.42 


s 


SW 5 


x h 


3/ 

- fa 


1/ 3/ 

/a - fa 


79.15 


79.14 ± 0.02 


-0.01 


0.29 ± 0.02 


2.84 ± 0.38 


0.54 


s 


SW 5 










79.14 ± 0.02 


-0.01 


0.29 Fix 


2.86 ± 0.21 


0.43 


s 


SW 5 










79.17 ± 0.07 


+0.02 


0.56 ± 0.07 


2.22 ± 0.80 


0.89 


s 


SW 4 










79.14 ± 0.04 


-0.01 


0.29 Fix 


1.47 ± 0.34 


0.51 


s 


SW 4 


fa 


3/ 
- fa 


9/ 7/ 
fa - fa 


65.21 


65.18 ± 0.03 


-0.03 


0.24 ± 0.03 


1.51 ± 0.45 


0.47 


s 


SW 3 










65.18 ± 0.03 


-0.03 


0.29 Fix 


1.64 ± 0.26 


0.30 


s 


SW 3 


3/ 
fa 


3/ 
- /2 


11/ 9/ 

fa - fa 


53.30 










0.73 


s 


SW 2 



Note to the table: A is an average wavelength for the various fine structure lines. 




200 400 600 800 1000 1200 1400 
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Fig. 7. Rotation diagram for the CO lines observed with 
the LWS (diamonds, this work) and for ground based data 
(filled circles), taken from Davis et al. (1999), Hogerheijde 
et al. (1999) and White et al. (1995). More than one LWS- 
value for the same upper energy E u /k refer to measure- 
ments with different detectors. The data have been (ad 
hoc) fit by two linear segments, the solutions of which are 
given in the figure (see also the text). 



cm 3 , which is in good agreement with the results by 
McMullin et al. (2000). 

In these models, the presence of a diffuse radiation 
field is introduced by the dust temperature T(j ust = 40 K, 
the wavelength of unit optical depth A r= i = 200 /im, 
the frequency dependence of the dust emissivity (3 = — 1 
and a geometrical covering factor of 0.5 (cf. Larsson 
et al. 2000). The (clearly detected) high- J lines are all 
only mildly sub-thermally excited (justifying a posteriori 
our initial assumption), but have substantial opacity, e.g. 



" : log JV co /dv = 17.3 and log(X c0 /dv/dr) = -7.4 




Rotational Quantum Number, J 

Fig. 8. The fluxes of rotational lines of CO, observed 
with the LWS, are compared to LVG model computations. 
Filled circles with error bars refer to LWS data, where 
more than one value for a given J correspond to different 
detectors. The upper limit is 3a. In addition, open sym- 
bols represent ground based data, where the different sizes 
refer to different telecope beams (see Fig.0). The physical 
parameters of the LVG model are indicated in the figure. 



7o(j=i4_i3) = 1-7. First at J u — 22 start the lines to be- 
come optically thin again (r (j = 22-2i) = 0.14). 

The principle parameter of the LVG model is related to 
the ratio of the column density to the line width, TV/ Aw. 
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For a given density of the collision partners, n(H 2 ), this 
ratio is given by 



Nco Y ,dv 
Av or 



where the right hand side is the 'LVG-parameter'. 
There, dv/dr is the (Doppler) velocity gradient in the 
gas and Xqo is the molecular abundance relative to H2. 
From Eq. (3), it is clear that LVG models are, in general, 
not unique, since an increase of the column density could 
have the same effect as a decrease of the line width. 

From the model fit, N(CO) = 1.5 x 10 18 cm~ 2 for 
the adopted Av — 7.5 km s -1 (FWHM of a Gaussian line 
shape; see Sect. 4.3). A circular source would have a diame- 
ter of about 5" (1500AU), a thickness of about 600 AU and 
an H 2 mass of about 0.01 M Q (for X CO = 10~ 4 ). Finally, 
the total CO cooling rate amounts to 3.6 x 10 _1 L©. 

The hot regions emitting in the H2 lines (Sect. 4.2.1) 
are not expected to contribute significantly to the CO 
emission observed with the LWS. We predict the strongest 
CO lines from this gas to be the ( J = 4 — 3) and the 
(J = 5 — 4) transitions, with "LWS" -fluxes from a 10" 
source of about 1 x 10 _14 erg cm~ 2 s^ 1 . 

4.2.3. LVG models: 13 CO, H 2 and OH 

The CO-model of the previous section can be used (by 
keeping dv/dr constant) to investigate whether it is ap- 
plicable also to other molecular species. Such a 'one-size- 
fits-air model would have the advantage of permitting the 
straightforward estimation of the relative abundance of 
these species (see Liseau et al. 1996b for an outline of 
this method). The reasonably satisfactory result of such 
computations for 13 CO, ortho-H 2 0, para-H 2 and OH is 
presented in Fig.^. 

The 13 CO spectrum has been computed under the as- 
sumption that 12 CO/ 13 CO is as low as 40 (Leung & Liszt 
1976). The data are clearly consistent with this value, but 
the S /N is insufficient to conclusively provide a better de- 
fined value. Since the 13 CO lines are all optically thin, the 
cooling in these lines (1.3 x 10~ 2 L Q ) is relatively more 
efficient than that in CO (by almost a factor of two). 

The H2O model is based on considering 45 levels for 
both ortho- and para-H 2 0, including 164 transitions each. 
The radiative rates are from Chandra et al. (1984) and the 
scaled collision rates from Green et al. (1993). The model 
fit of the observed spectrum requires an o/p = 3 for H2O 
and the derived H 2 0-abundance is X(H 2 0) — lx 10~ 5 . As 
expected, the excitation is sub-thermal and the lines are 
very optically thick (e.g., r 0(ol79 5Mm) = 433, r 0(pl0lMm) = 
218). Both the 380 GHz ortho-transition (4 M - 3 2 i) and 
the 183 GHz para-transition (3i3 — 2 2 o) are predicted to 
be masing (tq = — 1). The total cooling rate due to water 
vapour is L(H 2 0) = 2.1 x 10 _1 £o, i.e. at the 60% level 
compared to the CO cooling rate. 
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Fig. 9. The fit (smooth line) of the LVG model to the 
observed spectrum (histogram). All molecules, i.e. CO, 
13 CO, ortho-H 2 0, para-H 2 and OH, are assumed to 
share the same density and temperature, viz. n(H 2 ) = 
1.6 x 10 6 cm -3 , Tkin = 320 K; for further model details, see 
the text. Line identifications are as those given in Fig.||. 

For OH, the Einstein A values were computed from the 
data provided by D. Schwenke^, who also gives the energy 
levels. The collision rate coefficients for 50 transitions were 
obtained from Offer et al. (1994). Since rates are avail- 
able for collisions with both ortho-H 2 and para-H 2 , the 
total rate was based on H 2 -o/p = 1. As before, the excita- 
tion is sub-thermal and the lines are optically thick (e.g., 
T o(ii9^m) ~ 170 in each line of the doublet). This refers to 
the derived, relatively high, value of the OH-abundance of 
X(OH) = 2x 10~ 6 (OH/H 2 O=0.2). The OH lines cool the 
gas as efficiently as H 2 0, viz. i(OH) = 2.0 x 10" 1 L Q . The 
model is overpredicting the 119 /im line flux (whereas the 
113 (im H 2 line is underpredicted), perhaps indicating 
a distribution of temperatures (and densities). However, 
these lines fall in one of the least well performing LWS 
detectors (LW2) and instrumental effects cannot be ex- 
cluded. 

Anyway, we have so far considered only models of a 
homogeneous source at a single kinetic temperature in a 
plane-parallel geometry. The relaxation of these, likely un- 
realistic, assumptions is the topic of the next sections. 



Available at ittp://george.arc.nasa.gov:80/~dschwenke/ 



12 



Bengt Larsson et al.: The ISO-LWS map of the Serpens cloud core 

Table 6. Main input parameters of the dusty torus model 




Fig. 10. Model geometry of the dusty torus of SMM 1 
(see Sects. 4.2.4 and 4.2.5). Different shades of gray show 
schematically the density falling off outwards. The radius 
of the compact dense torus is ~100 AU, whereas the outer 
radius of the envelope is based on our maps. 
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Fig. 11. Comparison of the observed SED of SMM 1 and 
the model of the dusty torus. The individual fluxes (see 
Larsson et al. 2000 for details) are labelled by different 
symbols, to distinguish between beams of different sizes. 
The model assumes that we observe the torus at an an- 
gle of 31°, relative to its midplane. The effect of beam 
sizes is shown by the vertical lines and by the difference 
between the dotted and solid lines in the model SED at 
mid-IR wavelengths. Whereas only the lower points of the 
vertical lines are relevant, we have connected them to the 
adjacent continuum by straight lines, to better visualise 
the effect. The SED for the equivalent spherical envelope 
is also shown, to illustrate the influence of the bipolar out- 
flow cavities. 
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4.2.4. 2D radiative transfer model of the dusty torus 

In our previous paper, we presented a self consistent ra- 
diative transfer model for the SED of SMM 1 (Larsson et 
al. 2000). For a simplified analysis and for a direct com- 
parison with previous spherical models of the object, we 
adopted spherical geometry of the dusty envelope. The 
model provided a good fit to the observations longward of 
about 60 /im, but resulted in too low fluxes in the mid- 
IR. As already noted in that paper, the spherical symme- 
try may not be a very good assumption for SMM1, the 
source driving the bipolar outflow. In this paper, we per- 
formed detailed modelling of the dusty object using our 2D 
radiative transfer code (Men'shchikov & Henning 1997), 
which enabled us to quantitatively interpret existing dust 
continuum observations and to derive accurate physical 
parameters of SMM1. In the next section, the density 
and temperature structure of the model will be used in 
a Monte-Carlo calculation of the CO line radiation trans- 
fer in the envelope. Our approach and the model geometry 
are very similar to those for two other embedded proto- 
stars: HL Tau (Men'shchikov et al. 1999) and L1551 IRS 5 
(White et al. 2000); we refer to the papers for more de- 
tails on the general assumptions, computational aspects, 
and uncertainties of the modelling. 

The model assumes that SMM 1 consists of an axially- 
symmetric (quasi-toroidal), dense inner core surrounded 
by a similarly-shaped 'envelope' (Fig.[l(J). A biconical re- 
gion of much lower density with a full opening angle of 
100° is presumed to be excavated in the otherwise spher- 
ical envelope by the outflow from SMM 1. The structure, 
for brevity called 'torus', is viewed at an inclination of 30° 
with respect to the equatorial plane of the torus. Main in- 
put model parameters are summarised in Table|^. 

As very little is known about dust properties in SMM 1 , 
we adopted a dust model very similar to that applied by 
Men'shchikov & Henning (1999) for HL Tau and by White 
et al. (2000) for L1551 IRS 5. The dust population con- 
sists of 4 components: (1) large dust particles of unspeci- 
fied composition, with radii 100-6000 /im, (2) core-mantle 
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Projected baseline (k/E) Projected baseline fk/E.) 

Fig. 12. Comparison of the model visibilities at 0.8 mm, 1.4 mm, 2.7 mm, and 3.3 mm with available measurements of 
Brown et al. (2000) and Hogerheijde et al. (1999). The upper and lower curves in each panel show the visibilities for 
two directions in the plane of sky, parallel and orthogonal to the projected axis of the model. 
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Fig. 13. Density and temperature structure of the model torus of SMM 1. Left panel: The total densities in the midplane 
and in polar outflow regions, and dust densities of different dust grain components (for only their smallest sizes). Right 
panel: The temperature profiles correspond to the midplane of the torus; they were obtained self-consistently from the 
equation of radiative equilibrium. 



grains made of silicate cores, covered by dirty ice mantles, 
(3) amorphous carbon grains, and (4) magnesium-iron ox- 
ide grains. The latter 3 components of dust grains have 
the same radii of 0.08-1 /im. The dust-to-gas mass ratios 
of the components are 0.01, 0.0005, 0.0068, and 0.0005, 
respectively. The first component of very large grains is 
present only in the compact dense torus (7- <120AU), 
where all smaller grains are assumed to have grown into 
the large particles. Note that although unknown proper- 
ties of dust generally introduce a major uncertainty in the 



derived model parameters, extremely high optical depths 
in SMM 1 make the model results not very sensitive to 
the specific choice of the grain properties, except for the 
presence of very large grains in the dense central core. 

In the modelling of the dusty torus, we fitted all avail- 
able photometry of SMM1, paying special attention to 
the effect of different beam sizes. Important constraints 
for the density structure were provided by the available 
submm and mm interferometry of the object. The model 
SED, compared to the observations in Figs. Ill], fits almost 
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every single individual flux in the entire range from the 
mid-IR to mm wavelengths. Note that it would be wrong 
to fit the observed data with the total model fluxes, be- 
cause the angular size of SMM 1 is generally much larger 
than the photometric apertures. In fact, the model demon- 
strates that the effect of beam sizes on the fluxes may be 
as large as an order of magnitude. 

Comparison of the model visibilities to the interferom- 
etry data shown in Fig.|l^ demonstrates that the model 
is also consistent with the observed spatial distribution of 
intensity. The visibilities indicate that there is a dense core 
inside of a lower density envelope. The radial density and 
temperature profiles of the model, are shown in Fig. |l3|. 
The innermost dense core has ap cx r _1 density gradient 
in the model, whereas the outer parts of the lower-density 
envelope have a steeper density distribution (p cx r~ 15 ). 
The temperature distribution was obtained in iterations 
as a solution of the energy balance equation. 

4.2.5. 2D radiative transfer of the molecular emission 

We have used the density and temperature distributions 
of this dusty torus model in combination with a Monte 
Carlo scheme to compute the radiative transfer of the CO 
lines, and its isotopomers, through the source. 

Observations of the Serpens cloud core in the J = 
2—1 transitions of the CO-isotopomers C 18 and C 17 
are present in the archive of the James Clerk Maxwell 
Telescope (JCMT). These potentially optically thin lines 
could trace the embedded core SMM1. Our disk model 
reproduces the observed line intensities of these low- J iso- 
topomers fairly well (Fig.[l4|). There, the averaged back- 
ground emission of the surrounding gas has been sub- 
tracted, in order to reveal the line profiles of SMM 1 itself. 
From the figure it is evident that the C 18 and C 17 lines 
are optically thick out to a point, where the temperature 
falls below 15 K and where substantial condensation of 
the CO gas onto dust grains occurs. This CO freeze-out 
was treated following Sandford & Allamandola (1993 and 
references therein), where the ice-to-gas ratio r\ is propor- 
tional to the dust density n, and to functions of the gas 
and dust temperatures, T g and Td respectively, viz. 



r\ oc n 



T i/2 e i/T d 



(4) 



In the inner ('core') regions, line opacities are very high 
so that the photons are essentially trapped. Therefore, the 
excitation temperature of the molecules follows the kinetic 
temperature of the H2 gas. Beyond the core, i.e. in the 
'torus', r drops quickly below unity for the high-J lines, 
but there, the dust starts to contribute significantly to the 
(small) opacity (Fig.|l5|). 

The CO lines falling into the LWS range are all formed 
in the inner, hotter parts of the source (r <§; 10 3 AU, 
T ^S> 10 2 K). This small line forming region is insuffi- 
cient to produce the observed flux levels, i.e. the model 
undcrpredicts observed high-J line fluxes by more than 
two orders of magnitude. Irrespective of the geometry, we 
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Fig. 14. Upper panel: For background emission corrected 
line profiles of low- J CO isotopomers, viz. C 18 (2—1) and 
C 17 (2-1), toward SMM 1 are shown as histograms. The 
observations were retrieved from the JCMT archive. The 
results from 2D-Monte Carlo radiative transfer calcula- 
tions for the disk/torus model are shown by smooth lines. 
The shown integrated line intensity refers to the model, 
for which adopted abundances are 12 CO:C 18 0:C 17 = 
1 : 550 : 5 and freeze-out of the molecules is treated in 
accordance with the work by Sandford & Allamandola 
(1993). Lower panel: The total line centre optical depth 
in C 18 (2 - 1) and C 17 (2 - 1), along the line-of-sight 
toward the central region of the source, is shown by the 
solid line. Similarly, the dotted lines display the line opac- 
ities and the broken lines the dust opacities. 
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Fig. 15. Same as in Fig.|l4], but for three high-J CO tran- 
sitions, which fall in the ISO-LWS spectral band. 



can conclude quite generally that the CO lines observed 
with the LWS do not originate from the central regions 
of SMM 1, be it an accretion disk, be it infalling gas (we 
have also computed 'inside-out' collapse models). 

For the excitation of this gas we need to consider alter- 
native mechanisms and, since outflows are known to exist 
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in this region, shock heating of the gas offers a natural 
option. Our temperature determinations for the molecu- 
lar gas (Sect. 4.2.1) are also consistent with this idea. 

4.3. Shock heating 

From the discussion of the preceding sections we can con- 
clude that the heating of the gas is most likely achieved 
through shocks. These shocks are generated by flows 
within the LWS beam. Comparing the observed and pre- 
dicted molecular line emission with the J-shock models by 
Hollenbach et al. (1989) and Neufeld & Hollenbach (1994), 
we find that these models are in conflict with our obser- 
vations. 

In Fig.[L6|, we compare our observations of rotational 
lines of H2, CO, H2O and OH with predictions of the 
C-shock models by Kaufman & Neufeld (1996). The mod- 
els for logno = 5.5 (cm -3 ) and i> s ~15 - 20kms -1 are 
in reasonable agreement with the for extinction corrected 
(Ay =12 mag) observed values for H2 (Kaufman & Neufeld 
used o/p = 3) and for a flux from 6" x 6" (1 CVF- 
pixel). For CO, the model fits the observations for an 
adopted circular source of diameter 11". To achieve agree- 
ment for H 2 0, the model fluxes would need to be ad- 
justed downwards by a factor of 2.5, whereas an increase 
by more than one order of magnitude (a factor of 12) 
would be required for OH. Evidently, OH is largely un- 
derproduced by these models, a fact also pointed out by 
Wardle (1999). If on the other hand the Wardle model 
is essentially correct, this would suggest that the ionisa- 
tion rate in the Serpens cloud core is significantly higher 
(up to 10 -15 s -1 ) than on the average in dark clouds, 
£ = (10 -18 — 10 -17 )s -1 . High X-ray activity is known 
to be present within the Serpens cloud core (Smith et al. 
1999 and references therein). It is conceivable that such 
a high ionisation rate could also have considerable con- 
sequences for the cloud chemistry and its evolution. For 
instance, a relatively higher H^" abundance could be ex- 
pected, the effects of which (in addition to the enhanced 
abundance of OH) may in fact have already been observed 
(e.g., HCN/HNC~ 1; McMullin et al. 2000 and references 
therein). 



4. 4. Summarising discussion 

Based on their 0.8 mm JCMT-CSO interferometry, Brown 
et al. (2000) obtained estimates of the size, mass and av- 
erage (dust) temperature of the disk of SMM1. The es- 
timated mass is larger and the size of the disk is smaller 
by one order of magnitude than what is required to ac- 
count for the observed level of line emission (Sect. 4.2.2). 
Unless the disk (extended atmosphere?) is heated to very 
much higher temperatures (by an as yet to be identified 
mechanism) than the 60 K determined by Brown et al., 
we find it unlikely that the molecular line spectrum of 
SMM 1 is of circumstellar disk origin. Our own calcula- 
tions (Sect. 4.2.4) confirm this conclusion. 
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Fig. 16. Comparison of our molecular line observations 
with the predictions of theoretical models of C-shocks 
(Kaufman & Neufeld 1996). The pre-shock density is al- 
ways \ogn = 5.5 (cm- 3 ). Solid and dashed lines are for 
models with u s = 15kms~ 1 and 20kms , respectively. 
The H2 models refer to the CVF pixel size, whereas the 
other panels are for a circular source of diameter 11". 
The scaling factors, necessary to bring the observed and 
model data into agreement, are indicated in each frame 
by 'x factor'. 



It is intriguing that the luminosity of the spheri- 
cal model of SMM 1 (71 L®, Larsson et al. 2000) is 
close to the 'magic number' of the classical main ac- 
cretion phase of solar mass stars (Shu et al. 1987). At 
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the elevated cloud temperature of the Serpens cloud core 
(~ 40 K, White et al. 1995), the isothermal sound speed is 
0.4 km s -1 and, hence, the (time averaged, cf. Winkler & 
Newman 1980) mass accretion rate corresponds roughly to 
Mice = lCr 5 M yr -1 , yielding L acc = 70 L Q , where we 
have used the mass-radius relationship of Palla & Stahler 
(1990). In this scenario, the age of SMM 1 would be about 
10 5 yr or less, depending on the details of the acquired 
mass of the (presumably deuterium burning) central core. 
Regarding the data presented in this paper, we find it 
however difficult to reconcile this accretion shock model 
with our observations. As concluded in Sect. 4.2.4, the ex- 
citation of the observed lines requires significantly larger 
volumes at elevated densities and temperatures. 

The H 2 observations are partially resolved and there 
exists no ambiguity as to where, with respect to SMM1, 
the emission arises (cf. Fig.||). These lines trace a col- 
limated outflow toward the northwest of SMM1, which 
is also seen in ro-vibrationally excited H2 line emission 
(Eiroa & Casali 1989; Hodapp 1999). In the graphs of 
Fig.|l6|, we have assumed that also the LWS lines origi- 
nate essentially at the location of the H2 spots (i.e. we 
have artificially introduced another factor of two for the 
fluxes). However, the dereddened data of Eiroa & Casali 
(with the Ay-value determined in Sect. 4.2.1) could po- 
tentially present an additional difficulty for the C-shock 
model (Kaufman k Neufeld 1996). The estimated 1-0S(1) 
line intensity would in this case be larger by more than two 
orders of magnitude than that predicted by the model. We 
cannot exclude at present, however, the possibility that 
the 1-0S(1) emission observed by Eiroa & Casali (1989) 
is essentially unextinguished. Photometrically calibrated 
data at higher spatial resolution would be required to set- 
tle this issue. 

The mechanical energy input by the flow is £ m cch = 
0.5 /^ gas ?7iH no Vg x area, which for a pre-shock density of 
log no = 5.5, a shock velocity v s = (15 — 20)kms~ 1 , and a 
5" source size yields L mC ch = (6.4—15.1) x 10~ 2 L . From 
the Kaufman & Neufeld (1996) C-shock model, this gas 
is cooled by H 2 at a rate of (1.0 — 4.7) x 10~ 2 L Q . From 
the LWS data, we inferred the total cooling rate through 
the lines of CO, 13 CO, H 2 and OH of 78 x 10~ 2 L G 
(Sects. 4.2.2 and 4.2.3), corresponding to 0.5% to 1% of 
the total dust luminosity. This is larger by factors of 5 to 
12 and it is thus not excluded that the shocked regions ob- 
served in the H 2 lines and those giving rise to the FIR lines 
are not the same. We reached the same conclusion on the 
basis of our excitation and radiative transfer calculations. 

The observed and background-corrected [O 1] 63 fim 
emission toward SMM 1 suggests a contribution also by 
J-shocks within the LWS beam (Sect. 4.1.2). Intriguingly, 
the derived dimensions are practically identical to those 
determined for the LWS-molecular emission, albeit exist- 
ing J-shock models do not predict the relative intensities 
correctly. At present, we can merely conclude that shocks, 
in general, provide a plausible energy input mechanism, 
although the details of the shock type(s) are less clear. 
We propose that predominantly slow shock waves in the 



dense medium surrounding SMM 1 provide the heating of 
the molecules we have observed with ISO, whereas dy- 
namical collapse is not directly revealed by our data. 

5. Conclusions 

Based on spectrophotometric ISO imaging with the LWS 
and the CAM-CVF of significant parts of the active star 
forming Serpens cloud core our main conclusions can be 
summarised as follows: 

• We find the emission in the [0 1] 63 /im and 
[C 11] 157 /xm fine structure lines to be extended in our 
8' x 8' map. The absolute intensities and their ratios 
can be explained in terms of PDR models, where a UV 
field of Go = 15 ± 10 is falling onto the outer layers 
of the dark cloud, where densities are of the order of 
(10 4 - 10 5 )cm" 3 . 

• Also the emission in rotational lines of H 2 and high- 
J CO appears (slightly) extended, but we cannot ex- 
clude the possibility that it arises from point sources 
in the field, viz. the Class objects SMM9/S68 and 
SMM 4. The maximum of the emission is observed to- 
ward SMM 1, the dominating far infrared and submm 
source in the Serpens cloud core. 

• The spectrum of SMM 1 contains numerous lines 
of CO, H 2 and OH. These lines are generally 
subthermally excited and optically thick and trace 
regions of dimensions ~10 3 AU (V1500 x 600 AU), 
where temperatures are above 300 K and densities 
above 10 6 cm~ 3 . The derived abundances, relative to 
H 2 , are for CO, H 2 0, OH and 13 CO, respectively, 
X nw i = (1, 0.1, 0.02, > 0.025) x 10~ 4 . The ortho-to- 
para ratio for H 2 is consistent with the high temper- 
ature equilibrium value (ratio of the statistical weights 
of the nuclear spins), i.e. H 2 0-o/p = 3. 

• The relatively high OH abundance is indicative of an 
elevated level of ionising flux in the Serpens cloud core, 
causing the ionisation rate to be C ^ 10 _18 s _1 , i.e. 
significantly higher than the average rate prevailing in 
dark clouds. Strong and active X-ray sources, known 
to exist in the cloud, could be responsible for this ra- 
diation. 

• The observed SED of SMM 1 is consistent with a model 
of a dusty torus with an outer radius of 14 000AU 
(45"). The torus is heated by a central stellar source, 
with a (somewhat arbitrarily) adopted effective tem- 
perature of 5000 K. A luminosity of 140 Lq is required 
to explain the observations. The total mass of the 
toroidal core of SMM 1 is 33 M Q . The derived visual 
extinction through the torus exceeds 2000 mag and the 
torus is optically thick up to mm wavelengths. 

• 2D modelling of the radiative transfer through the cir- 
cumstellar torus of SMM 1 reveals that it is highly un- 
likely that the observed molecular emission arises in 
the torus. The same conclusion can be drawn for mod- 
els of dynamical collapse ('inside-out' infall). 
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• The H2 data have been obtained at significantly higher 
spatial resolution than that offered by the LWS. These 
CAM-CVF observations trace the regions of pure rota- 
tional H 2 line emission. The H 2 -maxinia are observed 
to be displaced from SMM 1 and are situated toward 
the northwest, along the jet of outflowing material 
from this Class source. The temperature of ~ 10 3 K 
of this H 2 gas is indicative of the heating by relatively 
slow shock waves. 

• The comparison of our molecular line data with shock 
models suggests that shock heating with v s ~ (15- 
20)kms _1 along the outflow of SMM 1 is the most 
likely mechanism of molecular excitation, although the 
details of these shocks are less clear. 
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Appendix A: Reduction of the LWS data 

The photometric calibration of the LWS is primarily based 
on observations and models of the planet Uranus. The S/N 
in these observations is rather modest, which will affect the 
relative spectral response function (RSRF) derived from these 
data. When calibrating the 'science data', the registered pho- 
tocurrent is divided by the RSRF, propagating any uncertainty 
in the RSRF, which will ultimately lead to errors in the derived 
flux density, F\. 

The derived photocurrent after standard processing, result- 
ing in an SPD file and where identified instrumental peculiari- 
ties have been removed, is shown for the LW 5 LWS detector in 



the upper panel of Fig. together with the relative spectral 



response function scaled to the same mean value. From the fig- 
ure, it is evident that many features seen in the photocurrent 
are due to the detector and grat ing response. 

The middle panel of Fig. A.l shows again the photocurrent 



and the RSRF, but now after the subtraction of a wide boxcar 
smoothing function (continuum subtraction). The upper curve 
shows the result of the division of the photocurrent by the 
response function. This allows us to gauge the effects on the 
spectra from the narrow features of the RSRF. In the long 
wavelength regime of the LWS, and in LW 5 in particular, there 
are spurious absorption features at positions, where the RSRF 
is steep. In regions, where the RSRF has steep gradients, can 
already very small wavelength errors create large features that 
could be mistaken for spectral lines. A wavelength mismatch 
between the photocurrent and the RSRF will also introduce an 
overall lowered S/N. 

In the lower panel of Fig. |A.l[ a relative shift by less than 
a quarter of an resolution element (the original sampling rate 
was at four times the spectral resolution) of the photocurrent 
resulted in the considerable reduction of an apparent broad 
'absorption' feature near 172 /im, an overall better S/N and 
therefore a better definition of the spectral lines. 

All individual spectra for all ten detectors have been care- 
fully monitored for obviously anomalous features which (most 
likely) were introduced by the RSRFs. This was done for the 
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Fig. A.l. Upper panel: The observed wavelength de- 
pendence of the photocurrent (thick line) for the detec- 
tor LW5 of the LWS. The thin line refers to the grating 
response function (RSRF) of the detector, scaled to the 
same mean value as the photocurrent. Middle panel: 
The lower curves are the same as in the upper panel, after 
subtracting boxcar fits. The division of the photocurrent 
by the RSRF results in the upper curve, which is the ob- 
served flux density (for clarity a median fit to the con- 
tinuum has been subtracted). Lower panel: The same 
as in the middle panel, but with a small wavelength shift 
(< 0.15^m) of the observed photocurrent, before division. 
Evidently, the 'noise' is much reduced: compare, e.g., the 
spurious 'absorption' feature near 172 fim. in the middle 
and lower panels. 



RSRFs of both OLP 8 and OLP 10 and any such spurious fea- 
tures were, of course, corrected for. This does not guarantee, 
however, that no such false spectral features do still exist in 
our data, as we were very restrictive in our application of any 
wavelength shifts. 
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Fig. B.l. Downward rate constants for collisions of CO 
with para-H2: Upper panel: as a function of J, where 
the solid lines are for T = 5, 10, 20, 40, 60, 100 and 250 K 
from Flower (2001), and the dashed lines for T = 
100, 250, 500, 750, 1000, 1500 and 2000 K from Schinke et 
al. (1985). Lower panel: as a function of Tkin, with the 
data by Schinke et al. adjusted to fit those by Flower at 
400 K, expanding the rates for the low temperatures to 
T = 2000 K (see the text). These rates are labelled 7 1 in 
the figure. 

Appendix B: Collision rates for CO 

We used rate coefficents for collisions of CO with H2 which are 
based on values found in the literature but which have been 
extended to rotational quantum numbers J u = 40, although 
extrapolations to higher J is not excluded. 

For J u < 29 and for low temperatures, (5 — 400) K, we 
used the recent rates of Flower (2001; ortho-H2-CO and para- 
H2-CO; downward rates). For the higher temperature range 
of (> 400 - 2000) K, the calculations by Schinke et al. (1985; 
para-H2-CO; upward rates) were used. The matching between 
these data sets is roughly acceptable, but there exist disagree- 
ments (Fig. B.l), which reflect the differences in assumptions 



and computational methods (see the discussion of resonances 
by Flower). 

In order to arrive at a consistent set of collision rate con- 
stants for the hole range of temperatures, the Schinke et al. 
data (correctly transformed to de-excitation rates; see also: 
Viscuso & Chernoff 1988) were laterally shifted to fit the 
Flower data at 400 K. A satisfactory matching was, however, 
not really possible for the lowest transi tion s connecting to the 
ground state, see the lower panel of Fig. B.l . In the figure, these 
expanded rates are labelled 7 . 

These 7 1 data span J-values up to 20 and temperatures 
between 400 and 2000 K. For the same temperature intervall, 
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Fig. B.2. Rate constants for CO as a function of J: 
Upper panel: the solid lines refer to the rescaled data 
of Fig. B.l, i.e. 7 1 , whereas the dashes are for T = 
500, 750, 1000, 1500 and 2000 K from McKee et al. 1982 
(collisions with He, scaled to H 2 , and up to J = 32). 
Lower panel: the fitted ratios of these rates ( J u < 20) as 
a function of J u and, as an example, for the five selected 
temperatures. 

McKee et al. (1982, MSWG) have published calculations (He- 
CO; downward rates) for J-values up to 32 (Fig. B.2 upper 
panel). These rates (xl.37) were divided into the Schinke et 
al. rates and fit by a p olyn omial to correct the shape of the 
McKee et al. data (Fig. B.3 upper panel), viz. 



ln7 I /ln7 MSW o =a + bj 20 < J < 32 
to obtain a new set of rate coefficients 



II _ a + bj 

7 ^MSWG 



(B.l) 



(B.2) 



For the expansion to J = 40, the 7 11 were fit to a second 
order polynomial, viz. 



In7 n = c + d.J + e,J 2 



(B.3) 



the result of which was used for the extrapolations of the para- 
H2 rates up to J = 40, viz. 



7 111 = exp (c + dJ + eJ 2 ) 



■J >32 



(B.4) 



This provides us finally with the full se t of collision rate coef- 
ficents for J u — > for J u up to 40 (Fig. B.3). 

So far, we have considered rates only connecting to the 
ground state, i.e. 7j_o- Flower (2001) did provide rates for 
collision transitions between all level's, but for higher J-values 
and/or higher temperatures we don't have that information. 
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Fig. B.3. Upper panels: Applying the curvature cor- 
rections to the data by McKee et al. results in the rates 
named 7 11 and shown by the solid lines. The extrapola- 
tions of 7 11 to high J-values are indicated by the dashed 
lines. Lower panel: The finally adopted collision rate co- 
efficients for CO, 7 111 , for the hole range in temperature, 
T = (5 - 2000) K, and energy, J = - 40. 

If the kinetic energy of the collision partners on the other 
hand is large compared to the rotational energy spacing of the 
CO molecule, the other rate cofficients can be obtained from 
(Goldflam et al. 1977; McKee et al. 1982) 

~<j-j> = (2 J' + 1) x 

V- <or , iA ( J L J ' V a2 (B.5) 
xE L (2L + l) ^ Q Q Q <,/7l-o 

where the Wigner 3-j symbol designates the Clebsch-Gordan 
coefficients, which were computed exactly (see, e.g., Zare 1986). 
Further, the factor Al,j 

= l + i ? L Z > where a = 0.13 Bo Z sca (£) * (B.6) 
6 + {aJ) z \T J 

given by McKee et al. (1982), should correct for the de- 
creased accuracy when the energy splittings become larger. 
-Bo = 1.9225 cm -1 is the rotational constant (Lovas et al. 1979), 
3A is a typical scattering length, /i = 3.5 amu is the re- 
duced mass of the colliding CO and He particles and T is the 
kinetic gas temperature. 

Finally, the rates for the inverse transitions were obtained 
from the condition of detailed balancing, viz. 

(2J' + 1) 7J ,_ J = (2J + 1) 7J _ J , exp(-huj_j,/kT) (B.7) 

For 13 CO, we used the same rate constants as for CO, with 
the proviso that the mass difference potentially introduces an 
additional error on the collision rates (3.5%). 
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